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We show that the Renyi entropies of single particle, critical wave functions for disordered systems 
contain information about the multifractal spectrum. It is shown for moments of the Renyi entropy, 
S n , where \n\ < 1, it is possible to extract universal information about the multifractility of such 
systems. This is shown through a generic calculation and then illustrated through two example 
models. We find good agreement between our analytic formula and numerical simulations of the 
two test models. Our formalism is easily extendable to generic non-interacting fermion models. It 
is also suggested that recent experimental advances in measuring the multifractal spectrum might 
allow some moments of the Renyi entropy to be measured. 
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I. INTRODUCTION 

There has been a recent influx of ideas and tools from 
quantum information theory that are being used to char- 
acterize many-body condensed matter systems. Part of 
this interest stems from the fact that some subtle proper- 
ties of the underlying phase of matter, otherwise difficult 
to calculate, can be found quite naturally by partition- 
ing the density matrix of a ground-state wave function 
and calculating measures of entanglement. For exam- 
ple, in one dimensional critical systems, it was shown 
that the central charge can be determined via the scal- 
ing of the entanglement entropy which takes the form, 
S ~ c/31ogL where L is the size of the partitioned re- 
gion and c is the central charge characterizing the critical 
point 1 . In higher dimensional critical systems general re- 
sults remain unknown, though interesting behavior has 
been found in a few isolated critical systems^ - — 

Entanglement finds a natural place in the study of dis- 
ordered systems where the concepts of entanglement and 
de-localization go hand in hand. For example, the en- 
tanglement entropy at infinite-randomness fixed points 
displays a similar log L scaling behavior in which, by 
analogy, it was proposed that these models may have 
very unconventional values for the central charge^. It was 
also conjectured that while the entanglement entropy has 
a similar behavior as the translationally invariant case, 
the Renyi entropies are very different, and in fact con- 
stant. Recent studies of the so-called single-particle en- 
tanglement entropy in disordered systems, in contrast to 
the many-body entanglement entropy described above, 
have found an interesting scaling behavior, S ~ ati logL, 
where L is the system size and ot\ is related to the mul- 
tifractal spectrum through, ot\ = d n T(n)\ n= i where t(ti) 
are the scaling exponents of the inverse participation 
ratios (IPRs) (defined below)2ii2. In this work we ex- 
tend this result and study the Renyi entropy for non- 
interacting disordered systems. Whereas references^ and 
[Tol only examined the n = 1 moment of the Renyi entropy 



S n , we show that the other moments yield more informa- 
tion about the multifractal spectrum. Indeed more than 
the derivative of the multifractal spectrum at n — 1 can 
be found by looking at other moments. Specifically, for 
small n the scaling behavior of the Renyi entropy and the 
IPRs are identical. Since the multifractal spectrum for 
the critical wavefunction is universal at Anderson transi- 
tion, the Renyi entropy for small n is also universa l 11 ' 12 . 
Thus we show that as in the case of critical one dimen- 
sional systems, universal information can indeed be found 
through entanglement. 

The multifractal spectrum (MFS) 11 ' 13 , like the central 
charge for critical 1+1 dimensional systems, plays an im- 
portant role in characterizing non-interacting disordered 



system: 



,14.15 



More concretely, the MFS is the set of scal- 



ing exponents of the IPR, defined as P n = \ijji{x)\ 2n 
which is indexed by a real number n. It has been observed 
that in many critical disordered systems the IPRs have 
a spectrum of scaling exponents defined by P n ~ L _T (") 
characterized by the spatial dimension d and function 
r(n) which is the MFS. Generically, the MFS can be 
written as r(n) = D n (n — 1) where D n is a function 
of n that depends on the system. For a conventional 
band insulator D n = while in a metal D n = d. At a 
critical point, the D n can have a non-trivial dependence 
on n*£. The MFS and its symmetry properties play a 
central role in restricting the possible types of critical 
field theories one can construc t 15 ' 17 " —. An understand- 
ing of the multifractal nature of critical wave functions 
is hence important for an understanding of the localiza- 
tion transition. Additionally, advances in imaging have 
also made it possible to begin to observe the multifrac- 
tal spectrum in a number of disordered systems at the 
Anderson transition 20 ^— . Our interest, however, in the 
connection between the MFS and the Renyi entropies, 
is chiefly theoretical. While for the sake of argument 
we consider models where the multifractal spectra are 
known, in general this set of exponents is difficult to com- 
pute. We demonstrate that a much simpler calculation 
gives an accurate approximation to an interesting region 



2 



of the MFS. By an explicit calculation and two numer- 
ical examples, we show that the information content of 
critical wavcfunctions of disordered free-fcrmion models 
contains imprinted signatures of their MFS. 



II. RENYI ENTROPIES 

The entanglement measures we focus on are the stan- 
dard von-Neumann entanglement entropy and a closely 
related quantity, the Renyi entropy. Given a wave func- 
tion \^), it is possible to construct the density matrix 
as, p = |\&)($|. The conventional entanglement entropy 
is calculated by partitioning the system into two spatial 
subrcgions A and B, where only one subregion is ob- 
served, say region A. One then traces out all of the possi- 
ble configurations of region B to yield the reduced density 
matrix pa = Ttb P- The question is then asked how much 
information does A know about region B and an obvi- 
ous candidate for such a measure is the von-Neumann 
entropy of region A, S = — Tr pA^ogp^A i-e. the entan- 
glement entropy. A closely related quantity is the set of 
Renyi entropies, 

S n = ^-logTr p n A . (1) 

parameterized by a real number n such that as n — > 1 
the Renyi entropy reduces to the von-Neumann entropy. 

We are focusing on localization in terms of single par- 
ticle physics and thus we will use a notion of single- 
particle entanglement for a choice of a single-particle 
eigenstate. One can define entanglement using a site 
occupation number basis in the second quantized Fock 
spaced. From now on, we work with a single-particle 
eigenstate ip(x) which can be written in the single par- 
ticle occupation basis as done by previous author o 9 ' 10 . 
Explicitly, the wave function can be written as 



\tp}= V>(r)|l)r ®r-7*r' |°>f 



(2) 
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where ^{r) is the normalized single particle wavefunction 
and \n) r is the state with n particles on site r. The 
wavefunction can then be written as, 



where \l) ri ^A are states where the single particle is lo- 
cated in region A at site r^. The matrix My generically 
has the form, 
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aj, bi denoting the amplitude of the wavefunction at site 
i in region A, B respectively. The Schmidt decomposi- 
tion of such a wavefunction is easily carried out. Using 



singular value decomposition, we can define a new basis 
u\-\])a and Vij\j)s where the matrices u,v are the uni- 
tary matrix of singular vectors for the matrix M. That 
is, they satisfy the relationship 



M, 



(5) 



where Xkk is the diagonal matrix containing the singular 



values of the matrix M,;, 
ements in the single particle problem 

A 2 = y/YT 



There are only two non-zero el- 
A i = VHi b l and 



af. The reduced density matrix pa is then 
easily written since the density matrix is now diagonal. 
One simply finds that 



Trpl = 



E< 



E 6 * 

Lies 



Pa+Pb 



(6) 



where Pa/b is the probability that a particle in the single- 
body wavefunction lies in region A or B respectively. 
Hence, we may write the Renyi entropy as 



1-71 



log(pS+PS), 



(7) 



which agrees with the expression for the entanglement 
entropy obtained by Ref. [lOl as n — >• 1 . 

Since we have restricted our attention to the one par- 
ticle sector of the density matrix and not the full many- 
body density matrix, the entanglement entropy should 
not grow large with the subsystem size as has been noted 
previousl y 10 i 24 . The expectation is that the same is true 
here. All the interesting behavior is in the scaling with 
the total system size. Indeed, by looking at the entan- 
glement of a single site, we can extract the dominant 
scaling behavior with the system size, L. Looking at the 
case where region A is a single site Xi, we simply have 
Pa(xi) = xf where d is the dimension. Clearly, any site 
Xi is not especially privileged and hence, one should look 
at the average over all possible sites, chosen with a uni- 
form probability distribution. Denote [-] Xi as this site 
average with respect to the uniform probability distribu- 
tion. That is, in d dimensions each site Xi is chosen with 
equal probability (a/L) d where a is the lattice constant. 
If we look at the site averaged Tr p n A we obtain, 



(3) [Tr^k= - 



((l-PA^r+piixi)) 



a\ d 
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(-i) m Ef™+E^ 



(-l) ro P m + P n 



n \ r(tn) / a \ T~(n) 

1 1 + (l 



(8) 

The quantity Ylx - PA.( x i) can ^ e identified with the in- 
verse participation ratio, P n which we recall scales as 
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Pn ~ (cl/L) t ^ 1 \ where the r(n) are the multifractal spec- 
trum. One finds that for \n\ < 1 the multifractal spec- 
trum is strictly less than zero and the site-averaged Trp^ 
can be approximated as, 



[Tr p n A ] Xt ~ 1 - n{a/L) d + (a/L) 



T{n)+d 



(9) 



By extension, the site averaged Renyi entropy can be 
approximated as, 



(1 - n)S n 



L 



-(n)+d 



\n\ < 1 (10) 



which is the main result of this work and clearly shows the 
relationship between the Renyi entropies and the multi- 
fractal spectrum. Grassberger— found a related connec- 
tion between Renyi entropies and the multifractal spec- 
trum of the probability measures of classical strange at- 
tractors. 

As n — > 1, as expected, our result reproduces the result 
of Jia et. ali£ for the von Neumann entropy, up to a sub- 
leading correction that vanishes as L — > oo. Indeed, there 
is always a region where r(n) is less than zero where our 
approximation is valid since the multifractal spectrum is 
a non-decreasing, convex function with t(0) = —d, the 
spatial dimension, and r(l) = 0. For moments \n\ > 1, 
such a simplification can not be achieved and according 
to Eq.(8), we need to sum up all the terms to get the 
Renyi entropies. The connection between the moments 
of the Renyi entropies and the multifractal spectra for 
these large moments is less clear. We will now consider 
two well-studied models where the multi-fractal spectrum 
is known and compare our approximate form with the 
known results. 
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FIG. 1. Comparison between the scaling behavior of (1 — 
n)S„+nL~ 2 ~ L 13 and the multifractal spectrum r 9c (n)+2 for 
the random-flux Dirac model for system sizes L=[100,1000]. 
In figures (a), the system is in the strong disorder regime 
while in figures (b),(c),(d) the system is in the weak disorder 
regime. In both, we find good agreement between the scaling 
of the Renyi entropy and multifractal spectrum for exponents 
\n\ < 1. 



Depending on the value of g, its been shown that the 
multifractal spectrum for a Dirac fermion in a random 
magnetic field has two regimes. In the weak disorder 
regime, defined by */ 2ir / g = q c > 1, the multifractal 
spectrum displays a discontinuous behavior as a function 
of n that was found in Ref. flTTl: 



III. RANDOM FLUX MODEL 

We first consider a system of Dirac fermions confined 
to 2d in the presence of a quenched random magnetic 
field normal to the plane with vanishing magnetic flux 
on average. Such a system can be thought of as the 
continuum limit of various tight binding Hamiltonians 
like the Chalker-Coddington network models. Gen- 
eral considerations have predicted that the random mag- 
netic field Dirac model has an exact zero-mode critical 
wave function which is multi-fractal. As has been noted 
by numerous authors, the most useful property of this 
model is that the wavefunction can be calculated ex- 
actly for any realization of the random magnetic field 
B(x) = V 2 <I>(x) and is given by ip(x) = Z~ x exp(— $(x)) 
where Z is the normalization factor. Here the ran- 
dom magnetic field B(x) is assumed to be Gaussian dis- 
tributed with variance (disorder strength) g, i.e. P[<&] oc 
exp[— J d 2 a;[V<I>] 2 }i2iiLL2!. To construct the reduced 
density matrix we simply note 



2(1 - sign(n)/(? c ) 2 n, q c < \n\ 
2(l-n/q 2 c )(n-l), \n\ < q c 



(12) 



In the strong disorder regime, where q c < 1, the similarity 
between the Renyi entropy and the multifractal spectrum 
can also be seen for values of the exponents given by 



-2(1 



H) 

n/q c ) 2 , 



q c < \n\ 
\n\ < q c 



(13) 



In both the weak and strong disorder regimes Eq. [TU] 
is applicable as long as n < 1. In Figure Q] we plot 
a few values of n and compare the scaling behavior of 
the Renyi entropy and the multifractal spectrum in the 
strong disorder and weak disorder regimes for system 
sizes L = [100, 1000]. Good agreement between the two 
quantities is found. 



IV. PERIODIC-RANDOM BANDED MODEL 



PA(Xi) 



-2$(xj) 



(11) To see that Eq. [TO] is quite general for any single par- 

ticle wave function, we explicitly check our result in a 
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FIG. 2. The multifractal spectrum and the scaling of the 
Renyi entropy moments as a function of n in the PRBM. The 
scaling behavior is shown for n = 0.2, 0.4, 0.6 and n — 0.8 
for system sizes L=[100, 1000]. In figures (a) and (b), the 
systems is in the weak multifractal regime with b — 4 and 
6=1 respectively. In figures (c) and (d), where b = 0.4 and 
b = 0.1 respectively, we can see that the agreement is good in 
the strong multifractal regime. In figure (d), only n — 0.2, 0.4 
and 0.8 are shown for the sake of clarity. 



much simpler model: the periodic random banded model 
(PRBM). The model is defined as the ensemble of ran- 
dom Hermitian L x L matrices where the entires Hij are 
independently distributed Gaussian variables with mean 
zero and a variance that falls off as 



m 3 \ 2 ) 



i 



i + (N-jW 



(14) 



At a = 1 the model undergoes an Anderson transition 
from the localized a > 1 to delocalized a < 1 phases, 
for all values of b. At this point, the model shows key 
features of the Anderson critical point, namely eigenfunc- 
tions possessing multifractal behavior. Again there is a 
single extended wave function on which we will focus. We 
expect and confirm that the previous analysis holds for 
moments less than n < 1. However, unlike the previous 



case, an analytic expression for the wave function is not 
known so we match our approximation to the numeri- 
cal calculation of IPR scaling. To eliminate the effect of 
boundaries, we will look at the periodic version of this 
model where 



1 sin 2 (7rr/A) 
P (7r/iV) 2 



(15) 



This can be interpreted as describing a ID model with 
long range hopping that falls off as l/r a — . From our 
previous results, it is expected that the combination 
(1 — n)S n + n-^ should scale with the system size to some 
power, (j-) P, where (3 is predicted to be r(n) + 1 (note 
the change from r(n) + 2 to r(n) + 1 when compared with 
the random flux case due to the change in spatial dimen- 
sion). Systems of size L=[100, 1000] were considered. As 
shown in Fig. f2J [3 matches very well with r(?i) + 1. In 
both the strong multifractal regime and weak multifrac- 
tal regimes there is a range of values for which the scaling 
of the Renyi entropy and the inverse participation ratio 
are nearly identical. 



V. CONCLUSIONS 

In conclusion, we showed that looking at the extended 
wave functions in the single particle occupation basis 
yields a general relationship between the Renyi entropies 
and multifractal spectrum. Our results can be applied 
to a wide-range of models, two of which were explicitly 
shown in this work. We hope that the Renyi entropy 
formula offers a simpler and more efficient way to calcu- 
late universal pieces of the multifractal spectrum. The 
connection also opens the door to measuring Renyi en- 
tropics experimentally through the measurement of the 
multifractal spectrum via imaging. Finally, it would be 
interesting to explore to what extent the Renyi entropy 
and the multifractal spectrum are related for many-body 
ground state wave functions. 
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